Numerical study of transport through a single impurity 
in a spinful Tomonaga-Luttinger liquid 
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The single impurity problem in a spinful Tomonaga-Luttinger liquid is studied numerically using 
path-integral Monte Carlo methods. The advantage of our approach is that the system allows 
for extensive analyses of charge and spin conductance in the non-perturbative regime. By closely 
examining the behavior of conductances at low temperatures, in the presence of a finite backward 
scattering barrier due to the impurity, we identified four distinct phases characterized by either 
perfect transmission or reflection of charge and spin channels. Our phase diagram for an intermediate 
scattering strength is consistent with the standard perturbative renormalization group (RG) analysis 
in the limit of weak and strong backward scattering, in the sense that all our phase boundaries 
interpolate the two limiting cases. Further investigations show, however, that precise location and 
form of our phase boundaries are not trivially explained by the standard RG analysis, e.g., some part 
of the phase diagram looks much similar to the weak backscattering limit, whereas some other part 
is clearly derived from the opposite limit. In order to give a more intuitive interpretation of such 
behaviors, we also reconsidered our impurity problem from the viewpoint of a quantum Brownian 
motion picture. 

PACS numbers: 71.10.Pm, 72.10.-d, 73.23.-b 



I. INTRODUCTION 

Low-energy excitations in a one-dimensional elec- 
tron system are described by the so-called Tomonaga- 
Luttinger liquid (TLL),^ which is characterized by 
power-low decay of various correlation functions and 
the spin-charge separation.^ Recent progress in nano- 
fabrication techniques has enabled us to realize quasi- 
one-dimensional quantum structures, in which TLL be- 
haviors have been experimentally observed in, e.g., 
fractional quantum Hall edges'^ and single-wall carbon 
nanotubes.^ One way to highlight such peculiar behav- 
iors of TLL is to introduce a single impurity, which dra- 
matically influences the transport of the system: in the 
low temperature limit, a conducting channel turns either 
perfectly transmitting or insulating. 

The impurity problem in a TLL is effectively described 
by bosonic fields at the impurity, which is equivalent to 
the problem of a quantum Brownian particle moving in 
a periodic potential. One of the powerful ways of treat- 
ing such a complicated quantum system is to numeri- 
cally simulate the system using the path-integral Monte 
Carlo (PIMC) method.^ Near the phase transition, how- 
ever, a simulation based on primitive local updates gen- 
erally fails, and one calls for a more effective update 
method optimized for the system. In recent years, the 
efficiency of the PIMC method has been remarkably im- 
proved by the extension of the Swendsen-Wang cluster 
algorithm^ to quantum systems with continuous degrees 
of freedom, and the algorithm has been first applied to 
the phase transition in a resistance-shunted Josephson 
junction system.^" 

The impurity problem in both spinless and spinful 
TLLs has been originally treated using perturbative 
renormalization group (RG) methods in the weak- and 



strong-impurity limits. Whether a conduction chan- 
nel becomes perfectly transmitting or insulating at low 
temperatures is determined by the relevance of the cor- 
responding backscattering or tunneling process. In con- 
trast to the spinless case, in which RG analyses in 
the weak- and strong-backscattering limits seem to be 
smoothly connected, the phase diagrams of the spinful 
case show clear inconsistency in the two opposite limits, 
i.e., mismatch of RG flow, suggesting the existence of an 
intermediate (unstable) fixed point. Besides, the phase 
boundaries between conducting and insulating phases for 
charge and spin are expected to shift continuously as a 
function of the backscattering strength. Although one 
can see the essence of critical phenomena in this spin- 
ful system using the standard perturbative RG approach, 
there is little information about the phase diagram for an 
impurity with intermediate strength. In this paper, we 
adopt the effective PIMC simulations mentioned above as 
a non-perturbative approach to study critical phenomena 
in the intermediate region of impurity strength. 

Our spinful impurity problem can be also understood 
as the physics of re- combination of the electronic charge 
and spin, which are generally separated and propagate 
with different velocities in a bulk TLL. One can introduce 
the scattering effect of an impurity, say, by annihilating 
one physical electron from the right-going mode, simul- 
taneously creating another in the left-going mode.^'^ We 
mean by a physical electron, an original electron com- 
posed of both charge and spin degrees of freedom. If this 
scattering potential is relevant and grows stronger, the 
charge and spin degrees of freedom become no longer in- 
dependent and their motion acquires some correlation. 
At low temperatures, the effect of an impurity becomes 
dominant. As a result the charge and spin tend to prop- 
agate almost together, realizing a situation which we call 
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spin-charge re- combination, unless the difference between 
their bulk velocities are too large. Of course, if that dif- 
ference is large enough, the charge and spin could re- 
main nearly independent, and the spin-charge separation 
preserves. We will argue, based on our numerical re- 
sults, whether or not the electronic charge and spin re- 
combine depends complicatedly on the competition be- 
tween the impurity strength and the difference in the 
original charge and spin velocities. 

This paper is organized as follows. In Sec. II we 
state our single impurity problem in a TLL — the spin- 
ful case with intermediate backward scattering strength. 
In Sec. Ill, we give details about the path-integral quan- 
tum Monte Carlo methods employed in this work. In Sec. 
IV, wc present our numerical results and the phase dia- 
gram deduced from our data, and then we discuss them 
in comparison with the known RG picture. Some further 
interpretation is also given in the context of quantum 
Brownian motion. Sec. V is devoted to the summary. 



II. STATEMENT OF THE PROBLEM 

Let us begin with introducing our model — the sin- 
gle impurity problem in a spinful Tomonaga-Luttinger 
liquid. Using the standard bosonization technique, we 
first formulate it in terms of two bosonic fields — one for 
charge, and the other for spin. Then, we rewrite it in a 
form suitable for numerical analyses. We also briefly re- 
view what is known about our model in the standard per- 
turbative RG picture. We end this section by addressing 
what we will attempt to uncover throughout this paper. 



order backscattering process is given by 

H,=VoJ2 ^i-^(O) + H.c. 

= ?; COS (/>p(0) cos 0^(0). (2) 

Here v is a. parameter proportional to the scattering 
strength Vq. Note that (/>y(0)/7r denotes the number of 
charges (spins) for v = p (cr) in the x > part of the 
system. Since the scattering term (2) influences only the 
fields at the origin, we can integrate out the other fields 
away from the barrier. If we write 4'v{t) = (j)i,{x = 0,t) 
in the imaginary-time formalism, the effective action 
takes the form 



S = So + Si, 



1/ UJ. 
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S\=V I 6.T cos (f)p{T) cos 4>cr{T) 

Jo 



(3) 

(4) 

(5) 



where (/>y(a;„) denotes the Fourier component of 0i/(t) 
and Un = 2im/(} is the Matsubara frequency, is 
known as the dissipative term in the Cakleira-Leggett 
model, and can be expressed by a form of long-range 
interactions in r direction as 



^0 



which is used when we apply the cluster algorithm to the 
PIMC simulation (see Sec. III). 



A. The single impurity problem in a TLL — the 
spinful case 

Low energy excitations of interacting one-dimensional 
electron system with spin are density fluctuations of 
charge and spin, labeled with subscripts p and a respec- 
tively, and the Hamiltonian is written as 

u—p^a " L \ / \ / _i 

where and 9^ are bosonic fields and the spatial deriva- 
tive of one field is the canonical momentum of the other. 
Ki, < 1 for repulsive interaction while > 1 for at- 
tractive, and Uv is the sound velocity of the density fluc- 
tuation. We now consider a single symmetric impurity 
with finite reflection localized at the origin, and introduce 
backward scattering of electrons by its barrier. Using the 
bosonized representation of a fermionic operator Vrs for 
an electron moving in the direction r = R or L with spin 
s =t or J,, the Hamiltonian corresponding to the lowest 



B. Consequences of the perturbative RG, and the 
three-dimensional RG phase diagram in the 
{Kp , , v)-space 

In order to allow for a comparison of our numerical re- 
sults with the known analytic viewpoints, here we briefly 
review the renormalization group (RG) picture presented 
in Refs. 6,7. The standard perturbative RG analyses 
can be performed either for an infinitesimal initial value 
of the scattering potential v in the original model (3) 
or in the opposite limit, i.e., for infinite backscatterings 
in the dual model of (3). They both give a phase di- 
agram in the {Kp, K(^)-plaTic characterized by four dif- 
ferent phases, which correspond to different transport 
behaviors of the system in the limit T ^ 0: (I) both 
charge and spin are insulating; (II) charge is conducting, 
while spin is insulating; (HI) charge is insulating, while 
spin is conducting; (IV) both charge and spin are con- 
ducting. Are the phase boundaries between such four 
different phases dependent on the initial values of v7 Ac- 
cording to Refs. 6,7, the obtained phase diagram (Fig. 1) 
in the above two limits have, as expected, similar configu- 
rations, but the phase boundaries are not located exactly 
at the same position in the {Kp, K„)-plarLe. 
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Let us now ask a question, what happens if we start 
from an intermediate value of the scattering potential 
u? For such a value of v, one can in principle consider 
a phase diagram, analogous to the above two limiting 
cases, i.e., probably with the same four distinct phases, 
but phase boundaries shifted from the two limiting cases. 
Since the bulk quantities Kp and K„ are invariant under 
the RG transformation (because the barrier is localized at 
the origin), we usually focus on a straight line connecting 
{Kp,Kc,0) and {Kp,Kc, oo) for a given set of Kp and 
K„, and examine how a scattering potential associated 
with a particular phase scales in the RG transformation. 
The results of Refs. 6,7 show that there exists a domain 
in the {Kp, K„)--pla.ne, in which this scattering potential 
is irrelevant in the limit of t; 0, whereas relevant in 
the opposite limit v ^ oo, indicating the existence of 
a non-trivial fixed point at an intermediate value of v 
(this fixed point is shown to be unstable). By performing 
the RG analyses one step further by considering higher 
order perturbations, one finds non-linear RG equations,^ 
suggesting a non-monotonous RG flow. 




FIG. 1: Three-dimensional phase diagram in the {Kp,KcT,v)- 
space. The phase boundaries axe analytically obtained in the 
following two limits: for a weak impurity {v ^ 0), the phase 
boundaries are three straight linos Kp + Ka = 2, Kp = 1/2, 
and Kc — 1/2; for a strong impurity (v oo), the phase 
boundaries are a hyperbola Kf^^ + K^^ = 2 and two straight 
lines Kp = 2 and K^ = 2. Three dotted lines represent those 
values of {Kp,Ka) at which phase boundaries in the v = 
and V = oo plane coincide. 



To summarize, the standard RG approach, with the 
help of duality transformation, not only reveals the RG 
flow in the limit of weak- and strong-backscattering barri- 
ers, but, by extending the perturbative analysis one step 
further, it also gives us some hints about how different 
tendencies of RG flow in the two limits evolve and eventu- 
ally merge in the region of intermediate coupling. On the 
other hand, there is little hope to obtain further informa- 
tion on the RG flow in the whole parameter space, most 
of which belongs to the so-called non-perturbative regime, 
by simply elaborating such an analytical approach. In 
this paper, we instead appeal to a numerical method, 
i.e., by performing a PIMC simulation for the effective 



action 5* given in (3), we study directly transport proper- 
ties of spinful TLL with an impurity of non-perturbative 
backscattering potential barrier. 



III. SIMULATION METHODS 

In this section, we illustrate our numerical simulations. 
In order to eliminate critical slowing down at low tem- 
perature and carry out efficient simulations of the paths 
4>p{t) and (f>„[T), we implement local updates in Fourier 
space and rejection- free global updates following Refs. 
10,12. Note that the single impurity problem in a spin- 
ful TLL is equivalent to the overdamped limit of the 
Josephson junction system discussed in Ref. 12. By dis- 
cretizing the imaginary time into N time steps, we define 
= MJP/N) (j = 0, 1, • • • , iV - 1). Then So and Si 
can be rewritten as 

N/2 

v=p,cr k=l 
N-1 

Si = vAt COS <j)pj COS (8) 

o = /^.iVV4A: fc=l,2,...,7V/2-l 
'^"^ - \ K^N k = N/2 ' 

where At = (3/N and ^^fe = (f>l,N-k = J2j <t>uj e^^^. In 

a local update for /c 7^ 0, a new value of (j)uk is randomly 
generated from a normal distribution with the variance 
f7,y/£^ in (9) by means of the Box-Miillcr method. This 
local update is accepted with a probability 

p = min{l,e-^'^i}, (10) 

where ASi is the variation of the potential term (8). For 

the fc = component (/),yo, a new value is generated from 
a uniform distribution ranged from — tt to tt, and the local 
update is accepted again with a probability (10). 

A global update scheme should be designed so that 
optimized paths for a given potential are efficiently gen- 
erated. In the case of the double cosine potential (2), 
an optimized path near the phase transition typically 
spends most of the time in potential minima, and also 
has some kink structures connecting adjacent potential 
minima, i.e., {4>p,(pa) = (^pTr, n^Tr) with integers rip and 
Tier such that Up + Ua = odd for t; > 0. In order to gener- 
ate such kinks, we apply the Swendsen-Wang algorithm^ 
to update of the continuous field variable following 
Ref. 10. To this end, we introduce a relative field variable 
ip^j = — (j)^"^"^™ as shown in Fig. 2, where the refer- 
ence is appropriately chosen as described below. 
If we regard the sign of the relative field s^j = <fvj/\Vvj\ 
as a spin variable, the dissipativc term in (6) can be rep- 
resented as a kind of one-dimensional long-range Ising 
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model 



So — Kvjj' SvjSuj' , 

j<j' 

, 1 



(11) 
(12) 



where each site labeled by j corresponds to each time 
step Tj = JAt, and has two spin variables Spj and Sg-j. 
Due to the finite bandwidth cutoff, we should represent 
a Fourier series and restrict the sum up to a 
cutoff frequency as 



2\ip^j\\ip^j, 



N/2 
k=-N/2+l 



rU-3')k 



(13) 



A cluster is built by connecting sites with the bond 
probability determined by the dissipative term Sq. In ad- 
dition, the cluster is flipped with no rejection, if ^j,™"'"'"'' 
is appropriately chosen so that the potential term S\ re- 
mains unchanged after the cluster is flipped. In Fig. 2, 
such rejection-free cluster updates implemented in this 
paper are illustrated in (pp-cp^ planes, where only the j- 
th path fragment in a cluster is shown. 

In the upper panel (a), two mirrors are located along 
^^mirror ^ (^^ + 1/2)7^ and <j)„'^'"°' = {71^ + 1/2)71, where 
Up and rio- are integers. Sites are connected with bond 
probability 



Pn' = 



c{0,l 



(14) 



and both the fields (ppj and tpa-j in the cluster are se- 
quentially reflected with respect to the two mirrors, i.e., 

{(ppj,ipgj) {-ippj,ipgj) {-ippj,-(paj). Note that 
the connected spins s^j and s^ji in one channel are not 
necessarily parallel, which is different from the original 
Swcndscn-Wang algorithm. After this double- fi,eld clus- 
ter update, kink structures connecting nearest-neighbor 
potential minima are inserted efficiently. 

Another cluster update is illustrated in the lower panel 
(b) in Fig. 2, where only a mirror for the charge degree 
of freedom is located at ^p™'""'' = np-K with an integer 
Ho- In this cluster is constructed with bond prob- 

ability 



Ppjj' 



max{0, 1 — e 



-2k„ 



(15) 



and the relative fields (ppj in the cluster are reflected with 
respect to the mirror, i.e., {ippj,(pa-j) {—ippj,(pa-j)- A 
similar cluster update can be performed also for the spin 
degree of freedom. These single-field cluster updates in- 
sert kink structures between next-nearest-neighbor po- 
tential minima. 

Using the PIMC method described above, we can 
efficiently simulate the impurity problem in a spinful 
TLL. In this paper, we observe zero-bias conductances 
of charge and spin channels to directly study the trans- 
port phenomena at low temperatures. In the linear re- 
sponse regime, a dc conductance at finite temperature is 
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FIG. 2: Cluster updates of field variables 4>pj and aj in a 



4>p-4>cr plane, where ip^j = (ji^j — <j)^™"''°'' . The empty circles 
represent the potential minima, and the dashed lines represent 
the reference fields <^p°""°' and <)}„'^'"°\ Only the field vari- 
ables at the j-th. time step in a cluster axe shown, (a) During 
one period of a double-field cluster update, a point {^pj,<fcTj) 
is subject to mirror reflection twice, i.e., once with respect 
to (l)pj = 4>p^'"°\ and subsequently to 4>„j = 0<,""'"°''. (The 
whole process is (ippj,ipaj) {-ifpj,ipaj) {-ippj, -ip^j).) 
(b) As for a charge-field cluster update, ip^j is subject to mir- 
ror reflection with respect to <j)p = (^p™"'"'"'', i.e., {ifipj, ippj) 
( — ippj,ippj). Similary, ip„j in a spin-field cluster is updated 
as {fpj,(paj) {(fipj, -(Paj). 



obtained from analytic continuation 



= lim G^{iiOn), 



(16) 



where the conductance at a Matsubara frequency can be 
calculated from a correlation function as 



Gu{iuJn) = f dr((/..(T)<^.(0))e--^ (17) 

Measuring the temperature dependence of the dc con- 
ductances (16) for different sets of Kp and near the 

phase transition, we determine the phase boundaries in 
the intermediate region of the impurity strength v. 
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FIG. 3: (color online) Symmetric coupling case : Kp ~ Ka- 
Charge conductance Gp(iuJn) (upper) and spin conductance 
Gaiiuon) (lower) for different values of Kp, plotted as a func- 
tion of tj„. K„ is fixed at K„ = 1.0 iV = 50, 100, 200, and 
400, and only the first ten points are shown for each Trotter 
number A'^. From top to bottom, the values of Kp are 1.2, 1.1, 
1.05, 1.025, 1.0, 0.975, 0.95, 0.9, and 0.8. Conductance curves 
(of both charge and spin) show an upward bend with decreas- 
ing uori when Kp > 1.025, whereas they are bent downwrad 
when Kp < 0.975. 



IV. RESULTS AND DISCUSSION 

This section is devoted to presenting our PIMC results 
and the phase diagram deduced from our data. We first 
present and analyze our conductance curves for several 
given points on the {Kp, ifCT)-plane, and determine what 
kind of phase those points belong to. We then discuss 
the whole phase diagram, and in particular the form and 
the position of our phase boundaries in comparison to 
their counterparts in the standard perturbative RG pic- 
ture (available only in the weak and strong backscatter- 
ing limits) . In order to uncover the nature of our phase 
boundaries, we also attempt to give further interpreta- 
tions to them in the context of quantum Brownian mo- 
tion. 



A. Charge and spin conductances and their "flow" 
at low temperatures 



FIG. 4: (color online) DC conductances deduced from Fig. 3. 
Charge (upper) and spin (lower) conductances for different 
values of Kp are plotted as a function of inverse tempera- 
ture N.Kcr is fixed at K^ — 1.0. Symmetric coupling case 
{Kp ~ Krj). As expected from Fig. 3, both charge and spin 
conductances increase with decreasing temperature (increas- 
ing A'') for Kp > 1.025, whereas they decrease for Kp < 0.975. 



We show in this subsection our simulation results for 
V = 4:,At = 0.25, and N = 50, 100, 200, and 400. Here 
iV is inversely proportional to temperature as = P/ At. 

Symmetric coupling case : Kp ~ K^ — Let us first in- 
vestigate a domain of {Kp, K^) for which the phase dia- 
grams in the limit of weak and strong scattering barriers 
are smoothly connected. This happens when two cou- 
pling constants are symmetric, or isotropic (see Fig. 1). 
In Fig. 3, we focus on the K^^ = 1 line (on which the spin 
part is SU(2) symmetric) and plot the first ten points of 
Gi,{iu)n) as a function of for different values of Kp near 
the phase boundary. For a given Kp, results for different 
N, i.e., for different temperatures are superposed to form 
a bundle of curves. For both Gp{iuJn) and Ga{ioj), one 
can see that the curves for Kp > 1.025 are bent upward 
with decreasing w„ (in the limit of cun — > 0), while the 
curves for Kp < 0.975 are bent downward. Note also that 
for a given Kp, the slope of different curves composing 
the same bundle always becomes steeper with decreasing 
temperature (increasing N). When the data shows such a 
monotonous dependence on temperature, one can deter- 
mine the phase boundaries by simply identifying a turn- 
ing point at which the bend of G^(iw„) changes from up- 
ward to downward with decreasing a;„. We will see later, 
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FIG. 5: (color online) Asymmetric coupling case : Kp <^ Ka- 
Charge conductance Gp{iuj„) (upper) and spin conductance 
Gcr{ii^n) (lower) for different values of Kp, plotted as a func- 
tion of uj„. Ka- is fixed at K^ = 1.8. N = 50, 100, 200, and 
400, and only the first ten points are shown for each Trotter 
number A'^. Charge and spin channels behave quite differently 
in this parameter regime. Upper : Gp{iuJn) is plotted for 
(from top to bottom) Kp = 0.600, 0.550, 0.525, 0.500, 0.475, 
0.450 and 0.400. With decreasing Un, the charge conduc- 
tance shows either a monotonous upward (first three curves) 
or downward (last three curves) bend. Such a behavior re- 
sembles the symmetric coupling case (see Fig. 3). Lower 
: Ga{iio„) is plotted for (from top to bottom) Kp — 0.500, 
0.450, 0.425, 0.400, 0.375, 0.350 and 0.300. The spin channel 
shows a non-monotonous behavior when Kp = 0.400, 0.375 
and 0.350. For details, see also Fig. 6. 



however, that the temperature dependence of G^{iu}n) 
curves can become non-monotonous in the presence of a 
non-trivial fixed point. We have actually determined our 
phase boundaries by tracing the temperature dependence 
of dc conductances obtained from (16) following Refs 
10,12. In Fig. 4, we plot the dc conductance of charge 
and spin as a function of the inverse temperature N , 
which are obtained by extrapolating the first five points 
on each curve in Fig. 3 to cj^ — > 0. With decreasing tem- 
perature, the conductances shows a monotonous increase 
(decrease) when Kp > 1.025 {Kp < 0.975). Here, the 
charge and spin channels show a simultaneous transition 
from conducting to insulating phase in consistent with 
the RG results. Recall that in the RG picture (see Fig. 
1) the IV and I phases touch at (Kp, K„) — (1, 1) both in 
the weak- and strong-backscattering regimes, which sug- 



FIG. 6: (color online) DC conductances deduced from Fig. 5. 
Charge (upper) and spin (lower) conductances for different 
values of Kp are plotted as a function of inverse temperature 
N . Kc is fixed at K^ = 1.8. Asymmetric coupling case {Kp <^ 
Ka). Upper : The temperature dependence for the charge 
channel shows monotonous increase {Kp > 0.5) or decrease 
{Kp < 0.5) with decreasing temperature. Lower : The spin 
conductance behaves non-monotonously for Kp — 0.35, 0.375, 
and 0.4, e.g., for Kp — 0.375, the conductance increases from 
N = 50 to A^ = 100, whereas it decreases from N = 200 
toward A = 400. 



gests that the phase boundary at that point is a straight 
line independent of v in the {Kp, K^, w)-space. 

Asymmetric coupling case : Kp <^ K„ — Let us turn 
to a parameter regime in which non-trivial RG flow is ex- 
pected for a finite backscattering strength. Such a behav- 
ior is actually expected whenever the phase boundaries 
in the weak- and strong-impurity limits are not identi- 
cal, but occurs typically when two coupling constants 
are highly asymmetric, or anisotropic : Kp <C K^. In 
Fig. 5, we plot Gp{iujn) and Ga-{ii^n) with K^^ fixed at 
Kcr ~ 1.8, and for various values of Kp. The correspond- 
ing temperature dependence of the dc conductances is 
shown in Fig. 6. A careful reader might immediately 
notice that the charge and spin channels behave differ- 
ently in these two figures. Of course, the origin of the 
difference is the anisotropy between Kp and K^-, but 
let us look into more carefully how they are different. 
In Fig. 5, conductance curves for the charge channel, 
i.e., Gp{iujn) with Kp = 0.600,0.550 and 0.525 are bent 
upward with decreasing a;„, while the same curves for 
Kp = 0.475, 0.450 and 0.400 are clearly bent downward. 
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The temperature dependence of the dc conductance (Fig. 
6) shows a monotonous behavior similar to the case of 
K„ = 1 (see Fig. 4). 

On the other hand, the temperature dependence of 
spin conductance Ga{ii^n) is more pecuhar: for example, 
if one focuses on the conductance curves for Kp = 0.35, 
their slopes are upward at high temperatures, e.g., be- 
tween N =50 and 100, whereas the same curves have an 
opposite slope at low temperatures, e.g., between N = 
200 and 400. Such non- monotonous behaviors might be 
more clearly seen, if we look into the dc conductance in 
Fig. 6 for Kp = 0.350, 0.375 and 0.400. 

Similar crossover behaviors are observed whenever an- 
alyzing the boundary between phases I and III, and are 
also reported in the Josephson junction system studied in 
Ref. 12. From the RG viewpoint,^ the unusual temper- 
ature dependence of the spin channel derives from non- 
monotonous flows of V near the intermediate unstable 
fixed point. Since the precise location of such a non- 
trivial fixed point is unknown, one cannot immediately 
conclude that the spin channel is in the conducting phase, 
even if the conductance, e.g. for Kp = 0.450 or 0.500, 
tends to increase monotonously toward low temperatures 
up to = 400. Possibly, it might turn insulating at a 
certain lower temperature which is numerically inacces- 
sible. Since it is difficult to locate true phase boundaries 
at T = in the presence of intermediate unstable fixed 
points, we instead identify the phase boundary at T / 
by observing the temperature dependence near the lowest 
temperature N ^ 400. 



B. The phase diagram in the (J<"p, ii'o.)-plane for an 
intermediate scattering strength 

Repeating the analyses outlined in the previous subsec- 
tion for different sets of Kp and K^, near the phase tran- 
sition, we determine the whole phase boundaries in the 
[Kp, Ka)-\i\a.Tie. In Fig. 7, we show our phase diagram 
for a finite impurity strength v — A obtained from the 
PIMC simulations at inverse temperature (3 = 400ZiT. 
Due to the symmetry of the action (3) in terms of p 
and (T, the following discussion also holds true when the 
charge and spin degrees of freedom are interchanged. In 
that case, the phases II and III are, of course, exchanged. 

In Fig. 7 one can see three different phases: (Phase I) 
neither charge nor spin is conducting, (Phase III) only 
spin is conducting, (Phase IV) both charge and spin 
are conducting. The phase boundaries are shown by 
solid lines. If one compares them with the phase bound- 
aries obtained by the renormalization group (RG) anal- 
yses in the weak and strong backscattering limits, one 
can verify that all the boundaries are indeed located be- 
tween the two limiting cases. However, the way they are 
shifted from either of the limits is not uniform, i.e., for 
nearly isotropic interactions Kp ~ ifo-, our phase bound- 
ary at an intermediate coupling is much closer to the 
weak backscattering (WBS) phase boundary, whereas for 
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FIG. 7: (color online) Phase diagram obtained from our PIMC 
simulations for the impurity strength u = 4. The solid lines 
stands for the boundaries between the I, III, and IV phases. In 
order to ease comaprison with the RG results, we superposed, 
as a guide to the eyes, the WBS results (dashed straight line) 
for the I-III and I-IV boundaries, and the SBS results (dotted 
hyperbolic curve) for the I-IV and III-IV boundaries, on top 
of our phase diagram. 



strongly anisotropic interactions Kp <^ K^, our phase 
boundary between I and III phases is much closer to 
the strong backscattering (SBS) phase boundary. Sim- 
ilarly, our III-IV boundary is much closer to the WBS 
boundary. Correspondingly, the tricritical point, i.e., the 
meeting point of I, III, and IV phases, lies between its 
counterparts in the WBS and SBS limits. 

For nearly isotropic interactions Kp ~ K^, the ob- 
tained critical line for a finite v between insulating (I) 
and conducting (IV) phases looks much similar to the 
WBS boundary: Kp + K^ = 2. This can be understood 
as a result of the large conductances of the charge and 
spin channels (see, e.g., Fig. 4). Then, we can judge that 
our scattering potential v is relatively small. For strongly 
anisotropic interactions Kp -C K^,, the III-IV boundary 
also resembles its WBS counterpart: Kp = 1/2, while 
the I-III boundary does not. Interestingly enough, for 
a broad range of Kp, say, < Kp < 0.3, the latter 
boundary is almost superposed on the SBS phase bound- 
ary: K^ — 2, which contradicts a naive expectation from 
the presumably small scattering potential. Thus increas- 
ing anisotropy, the position of our phase boundary shifts 
from that of WBS to SBS counterparts. 

In order to clarify such rather unexpected behavior of 
the I-III boundary, here we discuss how vertical and hor- 
izontal boundaries appear in the {Kp, K^)-plane. As an 
example of these boundaries, let us recall the RG phase 
diagram in the WBS and SBS limits in Fig. 1. In the 
WBS vertical boundary appears between III-IV 
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phases, where the interactions are strongly anisotropic 
{Kp 4C K^) and, for K„ > 2, the phase diagram is 
characterized only by the value of Kp. Duo to the 
strong attraction K^^ ^ 1, the spin channel transmits 
through the impurity so freely that the weak scattering 
potential hardly influences this spin channel. In that 
sense, the spin mode is irrelevant and the transport of 
the system depends only on the charge mode, which 
we call a one-field situation. In the SBS case, on the 
other hand, a horizontal boundary appears between I-III 
phases, where again the interactions are fully anisotropic 
and, for < 1/2, the phase diagram is characterized 
only by the value of K^. Due to the strong repulsion 
Kp ^ 1, the charge channel scarcely go over the impu- 
rity, and is almost extinct. Then again, we see another 
one-field situation. In both the WBS and SBS cases, 
the vertical or horizontal boundary appears as a result 
of the strong anisotropy in the interactions, rather than 
the extreme values of v. 

We can now interpret the unexpected change of the 
I-III phase boundary observed in Fig. 7. The horizontal 
region in the I-III boundary derives from the occurrence 
of a one-field situation where transport of the system 
is characterized only by the spin channel. Thus, after 
integrating out the extinct charge channel, one can argue 
that the effective action is given by 

^''T.^^\M^r.)\'+v I drcosMr), (18) 

which takes the same form as the action of a single impu- 
rity problem in a spinless TLL with interaction parame- 
ter K = Kc/2.^ Our phase diagram implies that the I-III 
boundary in the SBS limit: K„ = 2, partially preserves 
its position for a broad range of w > 0, leading to the ro- 
bustness of the boundary. The I-III boundary in Fig. 7 
also shows a small deviation from K^- = 2 with increasing 
the value of Kp. In this parameter region, the pinning of 
charge degree of freedom is no longer complete, and we 
believe that the crossover from the one-field model (18) 
to the original two-field model (3) occurs. We will give 
further discussion on this behavior in the context of quan- 
tum Brownian motion in the next subsection. Note that, 
due to the dualily of the impurity problem in a TLL,^''' 
the III- IV boundary will also show a similar crossover for 
a large value of v. 

It should also be added that a model qualitatively 
equivalent to the single impurity problem in a spinful 
TLL is realized in a completely different context. Werner 
et al. have performed the PIMC simulations of a sys- 
tem with two Josephson junctions, and shown a phase 
diagram similar to Fig. 7 consisting of three distinct 
phases. Although the kinetic term derived from charg- 
ing energy Ec, which is absent in our model, does not 
change the essential nature of the phase boundaries, it 
nevertheless influences transport phenomena at low tem- 
peratures. In Ref. 12, the authors seem to assume that 
the system undergoes transition to a one-field model like 
(18) as soon as a channel enters an insulating phase, and 
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FIG. 8: Minima and maxima of the two-dimensional potential 
V cos 4>p cos The empty (filled) circles represent the poten- 
tial minima (maxima). The solid (dotted) arrow stands for 
an example of tunneling of a massles particle between min- 
ima through a suddle point (over a maximum). For strongly 
anisotropic friction, e.g., Kf^^ ^ KcT^, the diagonal process 
(through a suddle point) is suppressed so that the relevant 
process becomes the vertical one (over a maximum). 



that the position of the tricritical point is independent of 
the Josephson coupling strength Ej, which corresponds 
to the backscattering strength v in our system. As we 
have seen above, however, the original two-field model 
slowly crossovers to a one-field one, and so the tricriti- 
cal point should in general depend on Ej, or v. If our 
conjecture holds true also for the two-Josephson-junction 
system, the discrepancy in the tricritical point discussed 
in Ref. 12 would be resolved. 



C. Interpretation in the context of quantum 
Brownian motion 

In the previous subsections, we have seen the sys- 
tem not only undergoes transition between insulation 
and conduction for charge and spin channels, but also 
crossovers from the original two-field (3) to a one-field 
model like (18). In order to interpret the crossover more 
clearly, let us reconsider our previous results from the 
viewpoint of quantum Brownian motion. As is clear from 
the Caldcira-Lcggctt form of the action (3), our spinful 
single barrier problem is equivalent to two-dimensional 
dynamics of a massless quantum Brownian particle in 
a periodic potential. In this picture, the bosonic fields 
{4>p, (f>a-) play the role of particle's coordinates. As shown 
in Fig. 8, the potential u cos 0p cos 0^ , for f > 0, has 
minima (maxima) at (n^Tr, na-ir) with integral Up and no- 
such that Up + Tier = odd (even). Each minimum corre- 
sponds to a certain ground state where integral numbers 
of electronic charges and spins exist in the a; > part 
of the system. The dissipation strengths in the (j)p {(j)^) 
direction is proportional to Kp^ (K^^). 

For nearly isotropic dissipations Kp^ ~ K^^, the par- 
ticle at low temperatures tunnels from one minimum to 
another, usually through a saddle point between them, 
and only occasionally over a maximum. Although charge 
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and spin arc generally separated and propagate with dif- done using the facilities of the Supercomputer Center, 
ferent velocities in a TLL, tunneling through a saddle Institute for Solid State Physics, University of Tokyo, 
point recovers the spin-charge combined nature of a phys- 
ical electron. On the other hand, for strongly anisotropic 
dissipations, e.g., Kp'^ :$> K^^, the position of the most 
relevant tunneling process could be taken by the other 
If the friction in the 4>p direction is large enough to sup- 
press completely the tunneling through a saddle point, 
the particle can only go over a potential maximum in the 
(j)a- direction toward another minimum. In this case, the 
system is dominated by the one-dimensional action (18) 
describing the horizontal part of I-III boundary in Fig. 
7. If the friction in the direction is so small that the 
massless particle can move freely in that direction, the 
system depends only on the (j)p coordinate and one gets 
another one-dimensional action analogous to (18). 

If we now return to the original TLL picture, the 
crossover from a two-field to a one-field model depends on 
the anisotropy of interactions and the impurity strength 
V in a, complicated way. Moreover, the phase diagram in 
Fig. 7 shows a crossover behavior in the Kp-K„ plane 
for an intermediate anisotropy, where the effective action 
can no longer be written in such a simple form as (18). It 
is worth noting that in the small- and large-barrier limits, 
crossovers to a one-field model occurs in so small a re- 
gion that one cannot observe them in the phase diagram 
in Fig. 1, while the PIMC simulation does demonstrate 
that they could actually appear in a broad range of v. 



V. SUMMARY 

In this paper, we have studied the single impurity 
problem in a spinful TLL using the path-integral Monte 
Carlo methods. Measuring the temperature dependence 
of the charge and spin conductances, we have obtained 
the phase diagram characterized by perfect conduction or 
insulation of the charge and spin channels, which is con- 
sistent with the renormalization group (RG) results in 
the weak- and strong- impurity limits. We have also ob- 
served non-monotonous temperature dependence of con- 
ductances, which qualitatively supports the non-linear 
flows near non-trivial unstable fixed points predicted in 
the RG picture. The phase diagram obtained from our 
simulations for an impurity with intermediate strength 
shows unexpected shift of a phase boundary for strongly 
anisotropic interactions. By mapping the impurity prob- 
lem to a quantum Brownian motion picture, we have pro- 
posed an intuitive interpretation of this behavior from the 
viewpoint of crossover to a one-field model. 
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